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ABSTRACT 

Context, To evaluate site quality and to develop multi-conjugative adaptive optics systems for future large solar telescopes, character- 
ed . ization of contributions to seeing from heights up to at least 12 km above the telescope is needed. 

Aims. We describe a method for evaluating contributions to seeing from different layers along the line-of-sight to the Sun. The method 
Pj^ ' is based on Shack Hartmann wavefront sensor data recorded over a large field-of-view with solar granulation and uses only measure- 

ments of differential image displacements from individual exposures, such that the measurements are not degraded by residual tip-tilt 
errors. 

Methods. The covariance of differential image displacements at variable field angles provides a natural extension of the work of 
Sarazin and Roddier to include measurements that are also sensitive to the height distribution of seeing. By extending the numerical 
calculations of Fried to include differential image displacements at distances much smaller and much larger than the subaperture 
diameter, the wavefront sensor data can be fitted to a well-defined model of seeing. The resulting least-squares fit problem can be 
solved with conventional methods. The method is tested with simple simulations and applied to wavefront data from the Swedish 1-m 
C-j ' Solar Telescope on La Palma, Spain. 

Results. We show that good inversions are possible with 9-10 layers, three of which are within the first 1.5 km, and a maximum 
I . distance of 16-30 km, but with poor height resolution in the range 10-30 km. 

Conclusions. We conclude that the proposed method allows good measurements when Fried's parameter r is larger than about 7.5 cm 
for the ground layer and that these measurements should provide valuable information for site selection and multi-conjugate devel- 
opment for the future European Solar Telescope. A major limitation is the large field of view presently used for wavefront sensing, 
, leading to uncomfortably large uncertainties in ro at 30 km distance. 
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1 . Introduction A unique advantage of solar observations is the availability 

, of fine structure for wavefront sensing nearly everywhere on the 

Presently, the most commonly used method for quantifying ^ surface Thifj allowfj , ^ tele s t0 te with 

astronomical seeing is the differential ima ge motion monitor ad iye icfj g ^ lock Qn ^ same M obseryed 

! <£^ M )' foposed for ESO site testing by |Sarazin & Roddier j ^ cameraf , A ticular blem of wavefront sensing 

' I t? ' ^lT n A°V!l e ZT^ C / 1C ^° nS SU f geStl °? with solar telescopes, however, is that low-contrast granulation 

O . oflFnsd l&m The DIMM built for ESO consists of a single stmcmres must ag and ^ necessita 6 tes the use 

<T. i 35 cm telescope with a pupil mask to allow images of bright stars of a faM , field . of . view (FOV) for wave front sensing. An 

> , to simultaneous y be monitored through two small subapertures j tant conseque nce of this is that contributions from higher 

• . separated spatially on a CCD by means of a beam splitter The , ^ a d oyer afl ^ ^ increases ^ hd ht This 

, advantage of this instrument is that differential image displace- d des ^ sensitivit t0 high . altitude see ing. 
£ . ments can be measured without impact from telescope tracking wkhin ^ framework of a desi stud for the European 

errors due to e.g. wind load. By comparing to theoretical cal- ^ Xele (EST) a program for studying two potential 

culations of Fried, this allows accurate seeing characterization ^ Qne Qn La palma and Qne Qn Tenerifei hag been initiated 

in terms of the so-called Fried parameter rn or an equivalent t,, , . . , . . , . , . ,. ... .. e ... 

, , F . _ " , , y 1 , . The goal is to characterize the height distribution of contnbu- 

seeing disk diameter related to r via 0.98A/r , where A is tiofls tQ geei ^ ^ twQ ^ and tQ define re irements for 

the wavelength In the following, we refer all values of r to a a multi . conjugate adaptive optics (MCA0) system for the EST 

wavelength of 500 nm. „ . , „ „ . in good seeing conditions. At present, very little work has been 

, T A " mstrument similar to te DIMM the S-DIMM published on characterization of high-altitude seeing based on 

r U !, & ff 6r S ^ ; u eC ^ rS l 00 , ? ^ ^c^'u 118 wavefront sensors operating on solar telescopes. Measurements 

the 4-m Advanced Technology Solar Telescope (ATST) but us- of scintillation of sunli ht with a linear of detectors have 

ing one-dimensional differential motion of the solar limb, mea- beefl shown tQ bg sensitiye tQ ^ hd ht distribution of seein 

sured through two 5 cm subapertures. contributions dSevkorall 1993r. Ib^cM1991l l2002h . Because of 

the integration of contributions to scintillation over the large 

Send offvrint requests to : l schanner@a S tro . su . se l solid an 8 le tended by the solar disk, an array of detectors 



2 



Scharmer & van Werkhoven: S-DIMM+ height characterization of day-time seeing using solar granulation 



with fairly large baseline is needed to achi eve sensitivity up to a 
height of only 500 meters (lBeckerslll993l) . Instrumentation that 
utilizes this technique, ref erred to as 'SHAdow BAnd Ranging' 
(SHABAR; Beckers 2001) was built to characterize n ear ground 
seeing in connection with site testing for the ATST dHill et al.l 
2006). A similar instrument, but with much longer baseline 
(3.2 m) than that used for the ATST, is under construction for 
EST site testing (Collados, private communication). 

The first published Shack-Hartmann (SH) based measure- 
ments of high-altitu de see ing with solar telescopes are those 
of Waldmann et al. (2008; 2010 in preparation), based on the 
covariance of local image displacements from different sub- 
apertures of a wide-field wavefront sensor (WFWFS). However, 
methods that are based on covariances (or correlations) of ab- 
solute image displacement s are sensiti ve to tracking erro rs and 
vibrations in the telescope (Fried 1975). Waldmann et al. (2008; 
2010, in preparation) therefore also implemented modeling of 
compensation of such image displacements with a tip-tilt mir- 
ror and adaptive optics (AO), assumed to correct Zernike aber- 
rations. A technique related to that proposed below is the 
SLODAR (Wilson 2002) for night-time measurements of seeing. 
The principle of this instrument is to use the cross-correlation of 
image motion measured for binary stars with a SH wavefront 
sensor through different subapertures. To eliminate the impact 
of telescope guiding errors, the average image motion of both 
stars measured with all subapertures is subtracted from the im- 
age motions of individual subapertures for each exposure sepa- 
rately. However, this procedure also removes the common atmo- 
spheric wavefront tilt, averaged over the entire telescope aper- 
ture for each of the two stars observed . This introduces an a niso- 
planatic component of the wavefront (iButterlev et alJl2006l) . As 
shown by the same authors, this can be compensated for in the 
analysis of the data, but at the price of making the analysis con- 
si derably more c omplex than in the originally proposed method 
bv lWilsonld2002l) . 

In the present paper, we propose a method that is also insen- 
sitive to image motion from the telescope or residual errors from 
a tip-tilt mirror but quite simple to model. The present method 
relies on measurements of the covariance of differential image 
displacements (between two subapertures) at different field an- 
gles. As we shall see, this approach can be considered as a nat- 
ural extension of the DIMM and S-DIMM approach to include 
measurements that are sensitive to the height variation of seeing 
and we therefore refer to the proposed instrument as S-DIMM+. 
We show that S-DIM M+ data can be inverted using calculations 
made bv lFriedl d 19751) to model contributions from different lay- 
ers along the line of sight (LOS). The paper is organized as fol- 
lows: In Sect. 2, we describe the assu mptions ma de, the method 
and extend the calculations made bv lFriedl (LL975). In Sect. 3, we 
describe the optical setup and the algorithms used to measure 
differential image displacements and to compensate for noise 
bias. In Sect. 4, we describe the results of the first data obtained 
with the S-DIMM+ installed at the Swedish 1-m Solar Telescope 
(SST) on La Palma and in Sect. 5 we give concluding remarks 
about the proposed method and its limitations. 



2. Method 

We make the following assumptions: 

1 . A measurement of image position is equivalent to a measure- 
ment of the slope of the wavefront over the subaperture. 

2. A measurement of image position averaged over the FOV 
used is equivalent to averaging the wavefront Zernike tip or 
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Fig. 1: Layout of the relation between the wavefront sensor geometry 
and the contributions to the two differential wavefront slope (image dis- 
placement) measurements from a height h, The separation between the 
two subapertures (indicated as heavy black lines) is s and the relative 
field angle between the two sub-fields measured is a. Note the diver- 
gence of the 'beams' with height due to the non-zero FOV used for 
these wavefront measurements. 



tilt over the solid angle corresponding to the subaperture and 
the FOV. 

3. The atmosphere is modeled as consisting of N discrete (thin) 
layers. 

4. The contributions from different layers are statistically inde- 
pendent of each other. 

5. The turbulent fluctuations in refractive index giving rise to 
the inferred wavefront slopes are statistically homogeneous 
and isotropic within a given layer. These fluctuations can 
be characterized in terms of a structure function based on 
Kolmogorov turbulence. 

6. Propagation and saturation effects are negligible. For a 
discussion of these effects on DIMM measurements, see 
iTokovinin & Kornilovl(l2007l) . 

The present method relies on using individual exposures to 
measure relative image displacements of the same target ob- 
served through different subapertures. These restrictions are es- 
sentially measuring differential image displacements from an 
individual exposed image, errors from telescope movements are 
eliminated; by using the same target observed through differ- 
ent subapertures we can use cross-correlation techniques with 
arbitrary solar fine structure to measure relative image displace- 
ments. We shall use the covariance of two such differential mea- 
surements, taken at different field angles, to achieve our goal of 
characterizing seeing. Furthermore, we shall restrict the present 
analysis to measuring only relative image displacements that 
are either along (longitudinal or x-component) or perpendicular 
(transverse or y-component) to the line connecting the centers 
of the two subapertures. Because of the assumed isotropy and 
homogeneity of the turbulent seeing fluctuations, the problem 
is then one-dimensional such that the measured relative image 
displacements can only depend on the distance between the two 
measurement points. Without loss of generality, we can assume 
that one of the subapertures is located at the origin, whereas the 
second subaperture is located at a distance s from that subaper- 
ture. Furthermore, we assume that the field angle of the first 



Scharmer & van Werkhoven: S-DIMM+ height characterization of day-time seeing using solar granulation 



3 



measurement is zero, while that of the second measurement is 
a and that this tilt is along the x-component of the image dis- 
placements. For the first measurement, the measured quantity 
6x\ then corresponds to added contributions from the N layers, 
located at heights h^. 

N 

^(i,0) = ^Wi)-x„(0)), (1) 

n=\ 

and that of the second measurement corresponds to 6x2 is then 

N 

6x2(s, a) = ^(x„(i + ah„) - x„(ah„)). (2) 

«=i 

Because of the assumed independence of contributions from dif- 
ferent layers, the covariance between 5x\ and 6x2 is given by 

{6x16x2) = ^{(x„(s) - x n (0)) (x„(s + ah„) - x n (ah„))), (3) 

n=l 

where (...) denotes averages over many exposed frames. 
Evaluating the four terms one by one and taking advantage of the 
assumed homogeneity of the statistical averages at each height 
h n , we can rewrite Eq. (O as a combination of three variances of 
differential image displacements 

N 

{6x16x2)) = ^{(x n (ah„ - s) - x„(0)) 2 )/2 

n=l 

+ {(x„(ah„ + s)- x„(0)f)/2 - {(x„{ah„) - x„(0)) 2 >. (4) 

This form of the equation clearly shows the connection to DIMM 
(ISarazin & Roddierl[l990l) and S-DIMM measurements. When 
a — 0, the last term disappears, the first two terms are equal 
and we have a conventional DIMM, except that the contribu- 
tions from the high layers are reduced by the averaging of image 
displ acements from the relativ ely large FOV (see below). 

ISarazin & Roddierl (Q990) gave an approximate equation 
for estimating the variance of differential image displacements 
recorded with two subapertures of diameter D, separated by a 
distance s: 

{(x(s) - x(0)) 2 ) = 0.358,l 2 r 5/3 Zr 1/3 (l - 0.541(s/£>r 1/3 ), (5) 

for longitudinal image displacements and 

{(y(s)-y(O)) 2 ) = 0.358^ 2 r 5/3 D- 1/3 (l -0.811( 5 /D)- 1/3 ), (6) 

for transverse image displacements. Here, D is the subaperture 
diame ter and rp isF ried's parameter. Adopting their notation and 
that of Fried! d 19751) . we rewrite these equations as 

{(x(s) - x(0)) 2 ) = 0.358,l 2 r 5/3 Zr 1/3 /(s/D, 0), (7) 
{(y(s) - y(0)) 2 ) = 0358A 2 r Q 5/3 D-" 3 I(s/D, n/2), (8) 

where it is assumed that the two subapertures are separated 
along the x-axis. The function / is normalized such that it ap- 
proaches unity when s approaches infinity and is symmetric 
/(-s/D,0) = I(s/D,Q) and I{-s/D,n/2) = I(s/D,n/2). With 

1 We will consistently refer to the heights h„ of seeing layers above 
the telescope as if observations were made with the Sun at zenith. For 
observations made at a zenith distance z, heights appearing in equations 
below need to be divided by cos z to correspond to the actual distance 
between the telescope and the seeing layer. 



an assumed zero inner scale and an infinite outer scale for the 
turbulence, the function / can only depend on separation s, di- 
vided by the diameter D of the averaging (round) aperture. The 
approximations given by Sarazin and Roddier are reasonably ac- 
curate only for s > D/2. However, to implement the present 
method we need more accurate estimates of / also when s/D ap- 
proaches zero and we need to take into account the averaging 
effect of using a large FOV for wave front sensin g (see below). 
We have extended the calculations of lFriedl d!975l his Eq. (32)) 
of /(s/D, 0) and I(s/D,n/2) to include also s/D < 1 and s/D 
as large as 50, corresponding to a distance of 5 m with 10 cm 
subapertures. The results of these calculations agree perfectly 
with those of Fried when s/D = 1 but systematically deviate 
with increasin g s/D and differ by 1 0% when s/D — 10. As sus- 
pected by Sarazin & Rod dier ( 1990), this suggests numerical in- 
accuracies in the calculation of Fried. We confirm this by noting 
that the calculations of Fried show a slight decrease of longitu- 
dinal differential image displacements when s/D > 7.5, which 
is clearly not physical. We obtain similar behavior in our cal- 
culations when the step size of one of the integrating variables 
(u in Fried's Eq. 32) is too large and have decreased that step 
size by approximately a factor of 80 compared to calculations 
essentially reproducing the results of Fried. With this improve- 
ment, the calculations show monotonous increase of I(s/D, 0) 
and I(s/D,n/2) with s as long as s/D < 300 which is a much 
larger range than needed for the present analysis. We therefore 
believe tha t our calculations must be more accurate than those of 
iFriedl (119751) . 

For wavefront sensing with an extended target, the averag- 
ing area corresponds to the sub-pupil area only close to the tele- 
scope. The effective averaging area expands from the pupil and 
up by an amount that increases with the FOV used for wavefront 
sensing. A reasonable estimate of the effective diameter D e ff can 
be obtained by calculating the convolution of the binary aperture 
and the FOV used for wavefront sensing projected at the height 
h. At small heights, the convolved function will be unity over an 
area corresponding to the sub-pupil and gradually fall off outside 
the sub-pupil. By calculating the area of the convolved binary 
pupil and FOV and equating that to 7rD 2 ff /4, an effective diam- 
eter D e tf can be defined. A good approximation is to set D e fj to 
the maximum of D and <ph, where <p is the (average) diameter of 
the FOV. At large heights, D e jf is much larger than the sub-pupil 
diameter D. At a distance of 30 km and with a FOV diameter of 
<p — 5.5 arcsec, we have <ph„ = 0.8 m, which is 8 times larger 
than the subaperture diameter D = 0.098 m used for record- 
ing data. The averaging effect of the large FOV used thus has 
strong effect on the determination of ro from differential image 
displacements measurements for the higher layers. In particular, 
we should not expect the method to work well, or even at all, 
when D e ft is larger than ro. This 'cone effect' is similar to that 
of measur ements of the integrated scintillation from the entire 
solar disk dBeckersll 1993b but of much reduced magnitude since 
the SH measurements use a typical FOV of 5-6 arcsec diameter, 
whereas the solar disk subtends a diameter 30 times larger. 

We conclude that we can use the theory developed bv IFriedl 
(fl97l for modeling differential image displacements measured 
with the proposed method by simply replacing s with s + ah 
and the pupil diameter D with an effective diameter D e ff. To 
make the averaging area roundish, we simply apply an approxi- 
mately round binary mask when measuring image displacements 
from the granulation images with cross-correlation techniques. 
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Combining Eqs. (0J, (|7]l and ©, we now obtain 

N 

(6xi6x 2 )) = 2^ c„F x (s, a, h n ), (9) 
«=i 

N 

{8yi6y 2 )) = Y u c n F y {s,a,h n ), (10) 

where 

F x (s,a,h n ) = I((ah n - s)/D eff ,0)/2 

+ I((ah n + s)/D eS , 0)12 - I(ahJD eS , 0), (11) 



F y (s, a, h n ) = I((ah„ - s)/D eB , n/2)/2 

+ l((ah n + s)/D eS , n/2)/2 - I(ah„/D eff ,n/2). (12) 

The coefficients c„ are given by 

c„ = 0358A 2 r (h n r 5/i D eW (hnr U \ (13) 

expressed in terms of Fried's parameter ro. We can also express 
the results in terms of the turbulent strength of each layer, C^dh, 

C\dh = 0.0599 cos(z)/l 2 r 5/3 , (14) 

where C\(h„) is the atmospheric structure constant and dh is the 
thickness of each (thin) layer. We obtain the relation 

c n = 5.98D e «(h„y 1/3 C?,dh/ cos(z). (15) 

Figure [2] shows the expected theoretical covariance functions 
given by Eqs. (fTTT) and (fT2l i for single seeing layers at differ- 
ent heights. We have made these calculations by assuming a 
5.5 arcsec diameter for the FOV used by wavefront sensing. 
We assumed ten subapertures of 9.8 cm diameter each across 
the 98 cm SST pupil diameter, giving a total of 85 well illumi- 
nated subapertures within the pupil. The maximum field angle of 
46.4 arcsec corresponds to the WFWFS built for the Swedish 1- 
m Solar Telescope. The calculations were made for heights h - 
0.0, 0.5, 1.5, 2.5, 3.5, 4.5, 6.0, 9.5, 16 and 30 km. From these 
covariance functions we conclude that the S-DIMM+ should be 
able to distinguish seeing at the pupil (h = 0) and at a height of 
about 500 m above the pupil and that the angular resolution is ad- 
equate to allow measurements up to a height of about 20-30 km. 
However, the similarity between the covariance functions at 16 
and 30 km height clearly demonstrates poor height resolution at 
these heights. This is a direct consequence of the large FOV used 
for wavefront sensing. As can be seen in the figure, the minimum 
height for which ro can be measured is set by the maximum field 
angle while the maximum height for which meaningful inver- 
sions can be made is set by the diameter of the FOV dWilsonl 
120021) . A particular feature in the covariance functions shown is 
the tilted dark line, marking minimum covariance. This corre- 
sponds to s = ah, i.e., where the beam defined by (s = 0,a) 
crosses the beam defined by [s, a = 0) such that a = sjh — see 
Fig. Q] This tilted dark line in principle allows direct measure- 
ment of the height of a single strong seeing layer, if that layer 
has a small vertical extent. 




a 



Fig. 2: Theoretical covariance functions, given by Eqs. jilt and d 1 2b . 
as function of separation s between the subapertures (increasing up- 
wards in each sub-panel) and field angle difference a (increasing from 
left to right). The left set of panels shows the covariance for longitudi- 
nal differential image displacements, the right panel that of transverse 
differential image displacements. The covariance functions have been 
calculated for heights 0.0, 0.5, 1.5, 2.5, 3.5, 4.5, 6.0, 9.5, 16 and 30 km 
and are displayed in order from bottom to top. The lowermost panel 
shows the number of measurements for each (s, a), applied as a weight 
W(s, a) in the least-squares fitting procedure. The lowermost row in all 
panels corresponds to i - = 0, for which the covariance functions are 
zero. 
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3. Implementation 

3.1. Wavefront sensor description 

The wide-field wavefront sensor (WFWFS) is mounted immedi- 
ately under the vacuum system of the SST. Light to the WFWFS 
is fed from a FOV adjacent to the science FOV and deflected hor- 
izontally such that the WFWFS beam does not pass the SST tip- 
tilt and AO system. The WFWFS optics consists of a field stop, 
a collimator lens and an array with 85 hexagonal microlenses 
within the 98-cm pupil diameter, the layout of which is shown 
in Fig. [3] The microlenses have an equivalent diameter close 
to 9.8 cm. We record data through a 10 nm FWHM 500 nm 
CWL interference filter using a Roper Scientific 4020 CCD with 
2048 x 2048 7.4 x 7.4yum pixels, operating at a frame rate of ap- 
proximately 9 Hz and with an exposure time of 3.0 msec. The 
image scale is 0.344 arcsec per pixel. 

3.2. Image shift measurements 

The lower panel in Fig. [3] shows a subimage with granulation, 
recorded with the WFWFS through one of its 85 subapertures. 
The subimage shown consists of 170x150 pixels, corresponding 
to 58.4x5 1 .6 arcsec. In the left- and right-hand parts of that panel 
are indicated a roundish mask, outlining two subfields at differ- 
ent field angles a, each with a diameter of approximately 16 pix- 
els, or <p — 5.5 arcsec. The granulation pattern within a mask of 
this and another subimage from another subaperture is used to 
measure differential image shifts between two subapertures, as 
follows: 

With the mask M(x, y) at the first subimage fixed, the mask 
at the second subimage is moved by (m,n) pixels in the (x,y) 
image plane, while accumulating the squared intensity difference 
between the two images, 

6I 2 (m,n) = ^}M(x,y)h(x,y)-M(x + m,y + n)h(x,y)} 2 , (16) 

x,y 

where l\{x, y) and h(x,y) correspond to the two granulation im- 
ages and M is unity within its 16 pixel diameter and zero out- 
side. This defines 6I 2 (m, n) at a number of discrete pixels sur- 
rounding the minimum corresponding to the best fit differential 
image shift. To improve the estimat ed image shift, t he quadratic 
interpolation algorithm proposed bv lYi et all (1 19921 Eq. (10)), is 
used. Tests with this image shift measurement algorithm indicate 
that an RMS error of about 0.02-0.05 arcsec is achievable when 
ro is in the range 7-20 cm (Lofdahl, in preparation). Following 
accepted nomenclature, we refer to the procedure outlined above 
as a "cross correlation" technique and 6I 2 (m, n) as "correlation 
function", although strictly speaking, we are not measuring im- 
age shifts with a cross correlation technique. 

3.3. Reference image frame selection and averaging 

The low contrast (typically 3% RMS in good seeing) granula- 
tion images observed through 9.8 cm subaperture and the rela- 
tively large (about 2 arcsec) granules, combined with the need 
to use a small FOV for wavefront sensing represents a signifi- 
cant challenge for wavefront sensing. Smearing of the granula- 
tion images by telescopic or atmospheric aberrations may cause 
large errors in the measured differential image shifts, or even 
complete failures. To reduce the risk of image shift measure- 
ments that fail completely, we do not perform image shift mea- 
surements directly between two arbitrary subimages. Instead, we 
use reference images, selected as the subimages with the highest 




Fig. 3: The upper figure shows the layout of the 85 fully illuminated 
hexagonal subapertures within the 98-cm SST aperture, indicated with 
the circle shown. Also indicated are two subapertures (s rl and s r2 ), cor- 
responding to selected reference subimages with high RMS contrast. 
Arrows point to two subapertures (s = and s) for which differential 
image shifts are measured with cross-correlation techniques. The mid 
panel shows the corresponding subimages with masks (not to scale) in- 
dicating two subfields at field angles separated by a. The lowermost 
panel shows one of the 170x150 pixel subimages with granulation struc- 
ture. The roundish masks outline the size and shape of the FOV actually 
used for differential image shift measurements. 



RMS contrast. In Fig. [3] top panel, is indicated two subapertures 
(at locations s r \ and s r 2) corresponding to such selected subim- 
ages and their relation to the two subapertures (s — and s) for 
which differential image shifts are measured. In the mid panel 
of Fig. [3] the corresponding subimages and the subfields at two 
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field angles, outlined with the roundish mask (M), are indicated 
schematically. The arrows in the top two panels of Fig.[3]indicate 
how image shifts are measured. We first measure the differential 
image shift between reference subfield 1 (at s r \) and the subfield 
at s=0 and also the image shift between reference subfield 1 and 
the subfield at s. We then calculate the difference between these 
two measured image shifts. We can express this operation as 

8x{s,a) = [x r \{s r i,a) - x(0,a)] cc - [x n (s r i,a) - x(s,a)] cc 

= x(s,a)-x(0,a), (17) 

where x r \(s r , a) is the position of the reference image mea- 
sured with a subaperture at an (arbitrary) location s r \ on the 
pupil and [...] cc indicates a measurement obtained with the cross- 
correlation technique outlined above. This cancels out the shift 
of the cross-correlation reference image, removes effects of any 
telescope movements (tip-tilt errors) and produces the relative 
image positions (5x, Sy) needed for the analysis. To further re- 
duce noise, we repeat this process 4 times, using the best 4 
subimages as references, and average the results of these 4 mea- 
surements. We emphasize that this corresponds to averaging 4 
measurements of the same quantity and that this is done sepa- 
rately for each pair of subapertures and subfields at each field 
angle a. 

3.4. Zero-point references and time averaging 

The differential image shift measurements obtained with cross- 
correlation techniques suffer from a lack of absolute zero-point 
reference. Such an absolute reference is in practice impossible 
to define precisely and furthermore not needed. Instead, we rely 
on the differential image shifts to average to zero over a suffi- 
ciently long time interval. The WFWFS data collection system 
is set to record bursts of 1000 frames, corresponding to approxi- 
mately 110 sec wavefront data. In the data reduction, we assume 
that the seeing induced differential image shift averaged over 
110 sec is zero and subtract the average shift measured from the 
1000 frames of the burst. This averaging and bias subtraction is 
done separately for each pair of subfields (a) within each pair of 
subimages at s and 5 = and for each of the 4 cross correlation 
reference images. This ensures that the measured covariances do 
not contain products of averages. Subtracting bias individually 
for data from the 4 cross correlation references is not needed 
for the covariance functions (we could just subtract the average 
shift for the image shifts averaged over the four cross correlation 
reference images), but is needed to correctly estimate the noise 
bias, as described in the following section. 

3.5. Noise bias estimation and compensation 

A fundamental limitation of this technique is the strong ground- 
layer seeing during day-time and the relatively weak contribu- 
tions from higher layers. A weak seeing layer with ro = 35 cm 
contributes (35/10) 5 ^ 3 = 8 times less variance than a seeing layer 
with ro = 10 cm. The effective wavefront sensing FOV diameter 
D e ff is 80 cm at 30 km distance (Sect. 2). According to Eq. ( IT31 ). 
the variance contributed from seeing at 30 km distance is thus 
reduced by approximately a factor (0.8/0. l) 1 ^ 3 = 2 by the av- 
eraging effect of the 5.5 arcsec FOV used to measure image 
displacements. Measurements of weak high-altitude seeing are 
therefore quite sensitive to random errors in the measured im- 
age positions, in the following referred to simply as "noise". In 
particular, this noise produces bias in the measured covariances 



that in turn leads to systematic errors, unless the noise bias is 
compensated for. 

To reduce noise, we repeat each differential image shift mea- 
surement 4 times with 4 different cross-correlation reference im- 
ages (Sect. 3.3). The corresponding image shifts are written for 
i =1^1 as 

6xi(s, a) = 6x(s, a) + «,(s, or)). (18) 

Here 6x(s, a) is the true shift and m(s, or) is the noise of mea- 
surement i. The four measurements are averaged to reduce noise 
before calculating the covariance matrix. In addition, these four 
measurements are used to estimate and compensate for covari- 
ance noise bias. The measured covariance C™ in the presence of 
noise is 

C™(s,a) = ([5x(s,0) + n(s,0)][5x(s,a) + n(s,a)]) 
= (5x(s, 0) Sx(s, a)) + (n(s, 0) n(s, a)) 

= C x (s,a) + 5C x (s,a). (19) 

where n{s, 0) is the noise averaged from the four measurements, 

n(s, a) = - 2^ ni(s, a), (20) 
i=i 

C x (s, a) is the covariance matrix in the absence of noise and the 
covariance noise bias is 

6C x (s, a) = (n(s, 0) n(s, a)). (21) 

To evaluate the expression for this noise bias, we assume that 
the noise of each measurement is uncorrelated with that of other 
measurements, (n,-(j, 0) nj(s, a)) = 0, unless i = j. We then ob- 
tain 

6C x (s, <*)= — ^ <«,(*, 0) m(s, a)) (22) 
i=i 

To estimate the noise bias directly from the data, we subtract the 
first measurement 6x\(s, a) from the corresponding average of 
all four measurements and calculate the covariance. We obtain 
the quantity 

Di(s, a) = ((m(sf, 0) - n(s, 0))(ni(s, a) - n(s, a))) 

13 1 4 3 1 4 \ 

= / [(-m(j, 0) - - 2 °>] [(4"^' <*) - 4 Yj ni(s > a)] ) ■ 

> i=2 i=2 ' 

(23) 

Evaluating this expression, we obtain 




(24) 



Forming similar expressions Dj(s, a), with i=2,3,4, for the re- 
maining three measurements and adding these, we obtain 

4 4 
2 D,(s, a) = || £<«,(*, 0) ruts, a)). (25) 
/=i 1=1 

Comparing to Eq. d22l i, this gives the desired estimate of the 
noise bias 6C x (s, a) as 

8C x {s,a) = — (26) 
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This noise bias is estimated at each s, a individually, for each 
data set separately and subtracted from the measured covariance 
matrix C™. 

We note that Eq. ( 1231 is insensitive to any random errors that 
are common to (the same for) all 4 measurements. Thus, the pro- 
posed noise estimation method is useful to indicate the magni- 
tude of random errors in the measured image shifts, in particular 
in bad seeing, but is not sufficiently accurate to provide robust 
estimates of the noise bias. Simulations are needed for proper 
understanding of noise propagation effects. 

3.6. Least-squares solution method 

To obtain the unknown coefficients c„ in Eqs. (|9]l and (TTOl i, we 
can solve a conventional linear least-squares fit problem by min- 
imizing the badness parameter L, given by 



( r 
L = 2j [C x (s, a) - 2^ c n F x (s, a, h„)] 2 

N 

+ [C y (s, a) - 2^ c„F y (s, a, h„)Y 



W 2 (s, a), (27) 



with respect to the coefficients c n . Here, the weight W(s, a) 
equals the number of independent measurements for each (s, a) 
(Fig. [2] lowermost panel). As is evident from the layout of the 
microlens array (Fig. 01, a large number of independent sam- 
ples is possible for small separations s, but only a few samples 
are possible for large s. Similarly, a large number of samples 
are possible for small, but not large, field angles a. The weights 
W(s, a) in Eq. ( TZTT i properly balance the least-squares fit to take 
into account the strongly varying number of measurements for 
different s, a and the corresponding noise. This is of particu- 
lar importance for (weak) high-altitude seeing, the signatures of 
which are at small angles a, where W is large, 

Minimizing L leads to a linear matrix equation Ac = b for 
c n . However, this permits solutions with negative values for c n , 
which is clearly not physical. In order to restrict the solutions to 
yield positive values for c„, we make the variable substitution 



c„ = exp(j„), 



(28) 



(Collados, private communication) and solve the corresponding 
non-linear least-squares fit problem with respect to the param- 
eters y„. Fits to data obtained so far indicates excellent conver- 
gence properties of the implemented non-linear method. 

3.7. Height grid optimization 

Good height grids can be found by calculating the inverse of 
the matrix A, corresponding to the linear solution for c n , and 
choosing a height grid that minimizes its noise sensitivity (sum 
of squared elements of the inverse of A). Such optimizations 
show that we should be able to determine contributions from the 
pupil plane plus about 8-9 layers above the pupil with the low- 
ermost layer above the pupil located at a height of 500 m. The 
maximum height can be 30 km with a FOV of 5.5 x 5.5 arcsec, 
however the height below that must be located around 16 km to 
not cause high noise amplification. The height resolution with 
which we can determine seeing contributions from layers above 
10 km with this large FOV is thus strongly limited. Only with 
a smaller FOV can the height resolution at large distances be 
improved. For the inversions discussed in this paper, we used 
the height grid defined by h= 0.0, 0.5, 1.5, 2.5, 3.5, 4.5, 6.0, 



30.0 km 
16.0 km 
9.5 km 
6.0 km 
4.5 km 
3.5 km 
2.5 km 
1.5 km 
0.5 km 
0.0 km 





10 20 
Height (km) 

Fig. 4: Response of inversion code to a single thin seeing layer of unit 
strength (c„ = 1) located at variable height (top panel). The vertical bars 
represent the c„ values at different heights (nodes), returned by the in- 
version code. The lower panel shows the returned integrated turbulence 
strength (2 c„) at each height. 



9.5, 16 and 30 km. We tested this configuration with input see- 
ing layers at heights in steps of 250 m from to 30 km height, 
then solved for contributions from the 10-layer height grid de- 
fined above. For input heights that matched one of the 10-layer 
heights, the inversion recovers the input height and contribution 
perfectly. When the input height is in between two of the heights 
in the inversion model, the inversion responds by distributing 
the correct c„'s between the two surrounding layers such that the 
relative distributions are in rough proportion to the difference 
between the true height and the two surrounding heights in the 
inversion model. By constraining the c„ values to be positive, 
negative overshoot in adjacent layers is prevented. This is illus- 
trated in Fig. |U which shows the response of the inversion to 
a thin seeing layer located at variable heighfl Due to the coarse 
height grid in the upper layers, the integrated turbulence strength 
is overestimated by up to 12%. A denser grid would reduce that 
overestimate but also increase noise amplification. We conclude 
that the method should work well with good input data. 



2 This plot is similar to that showing the heig ht response fo r the 
multi-aperture scintillation sensor (MASS) (Tokovinin et al. 2003|). 
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Fig. 5: The figure shows r for the ground-layer (h = 0) versus the 
integrated ro for the 9.5, 16 and 30 km layers and demonstrates that all 
inversions that return small values for r at the highest layers correspond 
to poor ground-layer seeing (small ro at h = 0). 



4. First Results 

We recorded 20 bursts of wavefront data, each consisting of 1000 
exposed frames, between 8:40 and 15:00 UTC on June 26, 2009. 
Each burst of 1000 frames took approximately 110 sec to record. 
After estimating and removing noise bias from each 1000 frame 
burst, the noise compensated data was analyzed in blocks of 
250 frames, corresponding to about 27 sec of seeing data. While 
recording this wavefront data (without AO), the SST was used to 
record science data with its AO system in closed loop. The initial 
seeing quality was considered good. Around 1 1 UTC, the see- 
ing started to deteriorate rapidly and science observations were 
stopped at 11:03 UTC. 

Due to problems with data acquisition, several data sets (16- 
19) contained corrupt images and were not processed. 

A fundamental aspect of the system is under what seeing 
conditions reliable inversions are possible. Several inversions re- 
turn very low values for ro in the uppermost layers. Figure [5] 
shows ro for the ground-layer (h = 0) versus the integrated ro 
for the 9.5-30 km layers. The dashed vertical line corresponds to 
ro = 7.5 cm for h — 0, the dashed horizontal line to ro = 25 cm at 
the 9.5-30 km layers. This figure shows that the smallest values 
for ro at the highest layers are all associated with poor ground- 
layer seeing. We conclude that inversions that return ro values 
at h = smaller than 7.5 cm, or about 75% of the subaperture 
diameter, should be rejected and excluded these data sets from 
further analysis. 

We now take a closer look at the data and the inversions. The 
covariance functions, defined in Eq. (0, were evaluated at steps 
of 5 pixels (1.72 arcsec) with the maximum field angle limited 
to 46.4 arcsec. Noise bias was subtracted from the data, as de- 
scribed in Sect. 3.5. The so-obtained covariance functions for 5 
data sets recorded between 08:40 and 11:01 UTC (zenith dis- 
tance in the range 30-61.0 deg) are shown in Fig. [6] together 
with the modeled covariance functions. The data and fits shown 
in this figure and in Fig.|7]are based on covariance functions cal- 
culated along rows of microlenses and corresponding sub-fields 
only (see Fig. [3]), in Fig.[8]are shown plots of observed and fitted 
data based on covariance functions calculated along both rows 
and columns of microlenses and sub-fields. The height grid used 
with these inversions has nodes at 0.0, 0.5, 1.5, 2.5, 3.5, 4.5, 6.0, 
9.5, 16 and 30 km. Comparing Fig.|6]to Fig. [2] the top two mea- 




Fig. 6: Covariance functions for 5 data sets processed (00:1, 03:0, 07:3, 
09:3 and 15:3). Within each of the panels is shown the observed (top) 
and modeled (bottom) covariance. Longitudinal covariances are shown 
in the left column, transverse covariances in the right column. The sepa- 
ration between the subapertures s increases upwards and the field angle 
a increases to the right in each sub-panel. Compare to the theoretical 
covariance functions for layers at different heights in Fig. [2] Data and 
fits shown are based on covariance functions measured along rows of 
microlenses and sub-fields only. 



sured covariance functions indicate clear signatures of a ground 
layer plus one dominating high-altitude layer (note the slanted 
dark line in these panels). In Fig. [7] are shown plots of these 
measured and modeled covariance functions as functions of s 
for field angles a of and 36.1 arcsec, where the second field 
angle is large enough to give very small influence from high- 
altitude seeing. It is evident that the fits are in general good, but 
that the measured transverse covariance is often too strong rela- 
tive to that of the longitudinal covariance. 

The jagged appearance of the plots in Fig.[8]are due to differ- 
ences in measured covariances from rows and columns. This ap- 
pears to indicate turbulence that is not isotropic above the ground 
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Fig. 7: Plots of measured (solid) and modeled (dashed) longitudinal 
(left column) and transverse (right column) covariance as function of 
pupil separations s at field angles a = arcsec (upper curves) and 
a = 36.1 arcsec (lower curves). Plots for data sets in Fig.[6](00:l, 03:0, 
07:3, 09:3 and 15:3, top to bottom) are shown, using only covariances 
calculated along rows. 



layer. We have made inversions with and without the column co- 
variances included and established that the overall conclusions 
drawn in this paper are robust, even though systematic differ- 
ences between the two types of inversions exist. 

4.1. Noise measurements 

In Fig. [9] (left panel) is shown the variation of the covariance 
(upper curve) and noise covariance (lower curve, multiplied by 
a factor 10) with or, at a separation s = 0.098 m for data set 3:0. 
Note in particular, that the noise covariance at s — is about 
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Fig. 8: Plots of measured (full) and modeled (dashed) longitudinal and 
transverse covariance as function of pupil separations s at field angles 
a = arcsec (upper curves) and a = 36.1 arcsec (lower curves). Plots 
for data sets 03:0 and 09:3 are shown (top to bottom), combining co- 
variances measured along both rows and columns. Compare to Fig. [7] 




rO al 16-30km (m) 



Fig. 9: Noise covariance and related data. The first plot shows the mea- 
sured covariance (upper) and 10 times the measured noise covariance 
(lower) at s = 0.098 m as function of field angle a for data set 03:0. 
This demonstrates the rapid drop (and reduced impact) of the covari- 
ance noise with field angle. The second plot shows the relation between 
the measured ground layer seeing (r at h = 0) and the noise in the mea- 
sured image positions (square root of the noise covariance). The vertical 
dashed line corresponds to r = 7.5 cm at h = 0, confirming our pre- 
viously determined threshold for useful data. The third plot shows that 
the noise in the measured positions increases with anisoplanatism, here 
quantified as the combined r for the two highest (16-30 km) layers 
used in the inversions. 



6% of the measured variance at a — 0, which is similar to the 
expected contribution from a seeing layer with ro = 35 cm at 
h = 30 km when ro = 10 cm at h = (see Sect. 3.5). Note also 
that the noise covariance drops off with angle in a way that is 
reminiscent of high-altitude seeing. However, the noise covari- 
ance drops by a factor of nearly 8 at a — 1.7 arcsec, such that the 
noise covariance is negligible already at this field angle separa- 
tion. This rapid decorrelation of the noise occurs on a scale that 
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500 m altitude shows a similar trend of degrading with time, but 
with turbulence strength that is typically 8 times weaker than 
those at h = 0. The seeing layer at 1 .5 km is even weaker than 
that at 0.5 km. 

The time evolution of seeing at the two highest layers are 
shown in the lower row of plots in Fig. [10] Except for a few 
data sets, the seeing contributions are consistently small from 
the 30 km layer. During this period, the high-altitude seeing 
primarily comes from the 16 km layer. The combined ro for 
the 9.5, 16 and 30 km layers averages at 37 cm. The accuracy 
of these estimates of high-altitude seeing is difficult to assess 
without independent verification. Images recorded with the 
SST (diameter 98 cm) generally show small-scale geometrical 
distortion but only minor small-scale differential blurring over 
the FOV, except when the Sun is at large zenith distance or when 
observations are made at short wavelengths around 400 nm. 
This is consistent with ro values of about one third of the SST 
diameter or larger, in agreement with the values found from this 
data. As regards the actual height of this seeing layer, upper-air 
sounding data above Tenerife at 0:00 and 12:00 UTC for this day 
(http : //weather . uwyo . edu/upperair/sounding .html 
Guimar- Tenerife) show a temperature rise above 17-18 km 
altitude, indicating the location of the tropopause, and enhanced 
wind speeds between roughly 10.5 and 15.5 km, peaking at 
14 km, altitude. If the latter layer is where the high-altitude 
seeing originates and if this is at the same height above La 
Palma as above Tenerife, then the height of that seeing layer 
should be about 11.5 km above the telescope. This corresponds 
to a distance of 23 km for the first data sets recorded at a zenith 
distance of 60 deg. This is midways between the two uppermost 
nodes in our inversion model. 

However, the seeing estimates at 30 km are very uncertain 
and our confidence in these estimates rely on the applicability 
of the weights Wis, a), defined in Sect. 3, and on the noise bias 
estimation and subtraction method outlined in Sect. 3.1. Setting 
all weights W(s, a) equal to unity reduces the average estimates 
of ro at 30 km by about a factor of two, but has only a minor 
effect on ro at all other heights. Similarly, reducing the weight 
W at a = by a factor of two also strongly increases the esti- 
mated turbulence strength at 30 km but has a small effect at other 
heights. These and other tests indicate uncomfortably large un- 
certainties of our estimates of seeing at 30 km distance, most 
likely primarily related to the large FOV used. 

An indication, based on SST observations, of seeing contri- 
butions from intermediate heights is the absence of noticeable 
large-scale variations of image quality over a science FOV of 
typically 1 arcmin. This suggests that the dominant seeing is 
close to the ground layer and at high altitude and that ro is sig- 
nificantly larger than 30 cm for intermediate layers. Our data are 
consistent with this. The solid curve in Fig.fTTIshows the turbu- 
lence strength C^dh as function of height, averaged over all data 
from 26 June 2009. Contributions to the seeing from heights in 
the range 1 .5-6 km are obviously small for this data set. A con- 
spicuous feature in Fig.QT]is the increased turbulence strength at 
3.5 km, suggesting a weak seeing layer. This feature may not be 
real, but an artifact of too dense grid points in this height range. 
We repeated the inversions after replacing the two nodes at 2.5 
and 3.5 km with a single node at 3 km, but leaving the remain- 
ing nodes unchanged. The result is indicated with a dashed line 
in Fig. QT]and shows a smoother variation of C^dh with height. 



Table 1: Signal to noise ratio for covariance functions measured at a = 
0, a = 1.72 arcsec and for de-correlated sub-fields (a > (p). 



is similar to the diameter of individual granules, rather than on a 
scale that corresponds to the diameter of the FOV. Table 2 sum- 
marizes the S/N, defined as the ratio of the measured covariance 
(corrected for noise bias) to the noise covariance, for a few se- 
lected data sets. In this Table, covariances used to calculate S/N 
have been averaged over all s in the range 0.098-0.784 m. The 
small S/N for data set 04:2 primarily comes from the excellent 
ground-layer seeing (ro = 42 cm) for this data set. Data set 20:0 
corresponds to one of the data sets rejected due to bad ground- 
layer seeing. Data set 22:1 corresponds to relatively poor seeing, 
ro = 7.7 cm for h — and ro = 6.9 cm integrated over the atmo- 
sphere, but evidently the wavefront data correspond to excellent 
S/N at all angles a. 

The mid panel of Fig. [9] shows the variation of the image 
position noise, calculated as the square root of the noise covari- 
ance at s = 0.098 m and a — 0, with ro at h = 0. Quite clearly, 
the image position noise increases dramatically when ro for the 
ground layer is less than about 7 cm. This renders wavefront 
sensor measurements essentially useless when ro < 7.5 cm, as 
concluded already from Fig. [5] Finally, the last panel in Fig. [9] 
shows a relation between the image position noise and ro at the 
highest layers, clearly demonstrating that anisoplanatism from 
high-altitude seeing causes problems with the large FOV used 
for wavefront sensing. When ro is larger than about 45 cm, the 
RMS noise is about 0.05 arcsec but when ro is 35 cm, the RMS 
noise is doubled. Since our FOV corresponds to a diameter of 
40 cm at 16 km, where most of the high-layer seeing origi- 
nates, our results confirm the expectation that the FOV should 
be smaller than ro at the height for which ro is determined, else 
wavefront sensor noise increases rapidly with decreasing ro. 

We finally caution that noise bias subtraction has a signifi- 
cant effect on measurements of seeing for the highest layers and 
that our confidence in these estimates relies on the method for 
noise compensation outlined in Sect. 3.5. It is therefore of con- 
siderable importance to study noise estimation by means of sim- 
ulations. 

4.2. Discussion of results 

Figures[10]and[TT]summarize some of the results. These two fig- 
ures are based on fits to measured covariances along both hori- 
zontal rows and vertical columns, corrected for noise bias (Sect 
3.5). 

The top row of Fig.[10]shows the variation of the turbulence 
strength C\dh with time for the three near-ground layers (h - 0, 
0.5 and 1 .5 km). The dominant layer is obviously at h — 0, and its 
variation with time clearly indicates gradually degrading seeing. 
The large scatter is real and illustrates the intermittent nature of 
the seeing, when averaged over relatively short time intervals 
(here, 27 sec), for these particular observations. The seeing at 
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Fig. 10: Summary of results. The variations of the turbulence strength C\dh with time are shown at different heights. From left to right: Turbulence 
strength at 0, 0.5 and 1.5 km height (top row) and at 16 and 30 km (bottom row). The lower-right plot shows C\ integrated over the 9.5-30 km 
layers. Note the different scales on the y-axis for these plots. 



4.3. Extension of the method 

In the present paper, we have calculated covariance functions 
from measured x, y positions along rows and columns of subim- 
ages and subfields. This was done in order to model the x,y po- 
sitions in terms of purely longitudinal and transverse image dis- 
placements. The advantage of this approach is simplicity as re- 
gards modeling but the disadvantage is that covariances can only 
be calculated from pairs of subimages and subfields that are on 
the same row or column of subapertures. A more appealing ap- 
proach is to process all data without these restrictions. This can 
be done by calculating Fried's function I(s/D e ff, 9) for arbitrary 
angles 9, defining the angle between a line connecting the two 
subaperture and that of the subfields measured. The measured 
x, y- positions should first be rotated onto the line connecting the 
two apertures at the pupil plane. This angle rotates with height 
h. A subfield FOV defined by a separation s at the pupil plane 
and field angles a and /3 are projected at (s + ah, fill) at a height 
h and tan# = f3h/(s + ah) which at large heights approaches 
/3/a. Thus virtually every subfield and every subaperture needs 
a unique basis function for every height layer in the model. By 
pre-calculating 7(s/D e ff, 9) and storing the result as a table, these 
basis functions can be calculated by interpolation. 

5. Conclusions 

The proposed method is based on measurements of differential 
measurements of seeing-induced image displacements, making 
it insensitive to telescope tracking errors, vibrations or resid- 
ual er r ors from a tip-tilt mirror. The numerical computations of 
iFriedl (119751) can be used to provide the theoretical covariance 
functions needed, requiring very small amounts of software de- 
velopment and avoiding the need for calculation of covariance 
functions via numerical turbulence simulations. The finite FOV 



used for wavefront sensing with solar granulation is accounted 
for in an approximate way by defining an effective subaperture 
diameter D e ff, increasing with height. 

In terms of data collection, the proposed method is iden- 
tical to the SLODAR method dWilsonl l2002h . also employing 
Shack-Hartmann wavefront sensor data. However, the SLODAR 
method uses averages of measured image shifts from all subaper- 
tures for each of the (two) stars observed to eliminate the effects 
of telescope guiding errors. Eliminating the anisoplanatism in- 
troduced _b^thj^_avera£ing_requires a fairly elaborate analysis of 
the data (Butte rlev et al.ll2006r> . In terms of analysis, the present 
method appears simpler. 

Based on simulations and data processed from a single day of 
observations, we conclude that the proposed method combined 
with wavefront data over about 5x5 arcsec subfields allows con- 
tributions to seeing from about 9-10 layers, stretching from the 
pupil up to 16-30 km distance from the telescope. At distances 
up to about 6 km, measurements with good S/N and a height 
resolution up to nearly 1 km appears possible. The 5.5 arcsec 
FOV used for wavefront sensing leads to poor sensitivity to high- 
altitude seeing and strongly reduced height resolution beyond 
10 km. At a distance of 30 km, the FOV used corresponds to 
averaging wavefront information over a diameter of 80 cm. Our 
estimates of r are very uncertain at this height and clearly the 
FOV needs to be reduced for seeing measurements at such large 
distances to be convincing. We also have established empirically 
that detecting seeing from high layers with the present system re- 
quires ro to be larger than approximately 7.5 cm for the ground 
layer and that the FOV should be such that the corresponding 
averaging area is smaller than ro at the high (10-30 km) layers. 

An important limitation is wavefront sensor noise, leading to 
bias in the measured covariances. By using image position mea- 
surements from 4 cross-correlation reference images, wavefront 
sensor noise is reduced and residual noise bias is estimated di- 
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Fig. 11: Summary of results. Shown is C\dh at different heights, av- 
eraged over all measurements from June 26, 2009 (solid curve). Also 
shown (dashed curve) is the average results for inversions made with 
the nodes at 2.5 and 3.5 km replaced with a single node at 3 km. Note 
that the varying zenith distance is not compensated for in the averages. 



rectly from the data and compensated for. However, the method 
used for estimation of noise bias relies on random errors from 
the 4 measurements to be independent and the method is fur- 
thermore "blind" to random errors that are the same for the 4 
measurements. Noise propagation clearly needs further investi- 
gation. 

The present method relies on wavefront tilts inferred from 
displacements of solar granulation images measured using 
cross-correlation techniques. The accuracy of these techniques 
for wavefront sensing is under investigation by means of sim- 
ulations (Lofdahl, in preparation). A weakness of such conven- 
tional techniques is that the image displacement is assumed to 
be constant within the FOV used for cross-correlations. This 
corresponds to modeling the wavefront as having pure tip and 
tilt without high-order curvature. At the intersection of adja- 
cent subfields, this corresponds to discontinuous gradients of the 
wavefront. Cross-correlations with overlapping sub-fields lead 
to multiple values of the wavefront gradient where subfields 
overlap. In addition to leading to inconsistencies and poor es- 
timates of wavefront gradients, the use of conventional cross- 
correlation techniques should also lead to noisier measurements 
when differential seeing within the FOV is strong. A more sat- 
isfactory approach may be to use a 2D Fourier expansion of the 
image distortions, 6x = dx(a,/3) and 6y = 6y(a,/3) over the en- 
tire FOV and to use that representation to calculate the local 
gradients of the wavefront. The highest order Fourier compo- 
nents would be limited by sampling and the FOV. This would be 
quite similar to fitting SH wavefront data to low-order Zernike or 
Karhunen-Loeve expansions, but over a rectangular FOV instead 
of over a round pupil. A simpler approach may be to use the 16 
pixel FOV cross-correlations as initial estimates and then refine 
the estimate with a 8-12 pixel FOV, restricting the corrections of 
the image position shifts to be small. 

We believe that seeing measurements from the ground layer 
may be possible with the proposed method and appropriate noise 
bias compensation even when ro is somewhat less than 7.5 cm 
but that little or no information can be obtained from the higher 
layers in such conditions. The major problem is the relatively 
poor daytime seeing, limiting the quality and number of seeing 



estimates for the highest layers. Possibly, the Moon can be used 
as widefield wavefront sensor target for verification of day-time 
estimates of high-altitude seeing at night, but with obvious lim- 
itations in image scale, exposure time and S/N requiring careful 
consideration. 
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